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Aggregation Driven by a Localized Source 

P. L. Krapivsky^ 

^Department of Physics, Boston University, Boston, MA 02215, USA 

We study aggregation driven by a localized source of monomers. The densities become stationary 
and have algebraic tails far away from the source. We show that in a model with mass-independent 
reaction rates and diffusion coefficients, the density of monomers decays as in d dimensions. 

The decay exponent has irrational values in physically relevant dimensions; /3(3) = (vTt-I- l)/2 and 
/?(2) = %/8- We also study Brownian coagulation with a localized source and establish the behavior 
of the total cluster density and the total number of of clusters in the system. The latter quantity 
exhibits a logarithmic growth with time. 

PACS numbers: 05.40.-a, 82.20.-w, 82.70.-y 


I. INTRODUCTION 


Aggregation underlies numerous phenomena ranging 
from polymerization [T] and aerosol behavior [5] to for¬ 
mation of clouds [3] and celestial matters where aggrega¬ 
tion plays a significant role in planet and star formation 
(and perhaps even leads to compositeness of dark matter 
a)- Earlier work on aggregation is described in Els]; 
the following work is summarized in a number of books 
and reviews [MU- The framework for studying dilute 
well-mixed aggregating systems was developed by Smolu- 
chowski who wrote an infinite set of rate equations for 
the concentration Ck{t) of clusters of “mass” k: 


dck 

dt 


2 X Cfc ^ Kk,jCj 

i+j=k j>l 


( 1 ) 


These equations reflect that in applications each cluster 
is often composed of an integer number of identical el¬ 
emental building blocks, ‘monomers’; then the number 
of constituting monomers can be called the mass of the 
cluster. Another tacit assumption is that the mass is 
the only relevant characteristic, so e.g. the shape of the 
cluster does not matter. This is certainly a simplifica¬ 
tion, although in some applications (e.g., when clusters 
are spherical or have some other fixed shape) it suffices 
to classify clusters only by their masses. 

There are only quadratic couplings between the den¬ 
sities in Eqs. Q since only binary collisions are taken 
into account. This is reasonable in dilute systems where 
binary interactions usually dominate |12| . The rate at 
which clusters of mass i and j merge is denoted by Kij . 
Mathematically, the reaction rates Kij form an infinite 
symmetric matrix, Ki j = Kjy, the pre-factor ^ in front 
of the gain terms on the right-hand side of Q is needed 
to avoid double counting. 

In Brownian coagulation 0181 n], the reaction rate 
admits a simple expression, Ki^j = ATT{Di + Dj){Ri-\-Rj), 
through radii and diffusion coefficients of the merging 
droplets. In the regime when droplets are sufficiently 
small, the Stokes-Einstein relation implies Di ^ l/i?i, 
and the Brownian coagulation kernel becomes 

= ( 2 ) 


Here we have dropped numerical factors and taken into 
account that Ri ~ for spherical droplet in three 
dimensions. Equations 0 with reaction kernel 0. are 
analytically intractable and Smoluchowski suggested to 
replace the Brownian coagulation kernel by the constant 
kernel, Kij = const. Equations 0 with this kernel are 
solvable. Apart from being the simplest possible kernel, 
the mass-independent kernel has the same homogeneity 
as the Brownian coagulation kernel ([^, namely in both 
cases K\iXj = Ki,j- Chief qualitative behaviors are sim¬ 
ilar in both cases, e.g., the typical cluster mass grows 
linearly with time: (fc) ^ t. (More subtle behaviors, e.g., 
the shape of the scaled mass distribution, differ pTHITT] .! 

Aggregation with homogeneous source has also at¬ 
tracted considerable attention, see [inillllllHIS]- The 
details of the source play rather limited role, e.g., they do 
not affect the emerging large mass behavior if only clus¬ 
ters of small mass are injected. It is customary to assume 
that only monomers are injected. The governing equa¬ 
tions then differ from Q by the term JSk,i, with J being 
the strength of the monomer flux. Since the source is as¬ 
sumed to be homogeneous, the governing equations are 
still ordinary differential equations. Aggregating systems 
with homogeneous source often reach a non-equilibrium 
steady state [ini HB nnn]; this makes such driven sys¬ 
tems more tractable than un-driven systems. 

Here we study aggregation processes with a spatially 
localized source. The densities Cfe(r,t) vary according to 
a system of partial differential equations (PDEs) 

DfPlJ Ck ~ ^ ^k ^ ^ d^k,j^j 

i+j=k j>l 

J6k,iS{r) (3) 

The first term on the right-hand side describes mixing 
due to diffusion, the next two terms account for aggre¬ 
gation, and the last term represents the monomer source 
at the origin. 

At first sight, the system (|^ of infinitely many coupled 
non-linear PDEs is totally intractable. We are chiefly in¬ 
terested, however, in the large time behavior when cluster 
densities may become stationary. Whenever this hap¬ 
pens, Eqs. 0 reduce to ordinary differential equations 
(ODEs). These ODEs are still generally intractable, so 



2 


one can try to get insight from the simplest model with 
mass-independent reaction rates supplemented by mass- 
independent diffusion coefficients. 

The remainder of this paper is organized as follows. 
In Sec. we analyze the model in which reaction rates 
and diffusion coefficients are constant. We begin with 
the total cluster density (asymptotic behaviors of the to¬ 
tal cluster density far away from the source were origi¬ 
nally studied in Ref. [3D]). We then derive more subtle 
decay laws for cluster densities. The results of Sec. [n| 
show that the basic quantities decay algebraically with 
distance from the source, and suggest to seek similar so¬ 
lutions for models where reaction rates and diffusion co¬ 
efficients are homogeneous functions of mass. In Sec. |III| 
we apply this strategy to Brownian coagulation with lo¬ 
calized input. Finally, in Sec. IV we give a summary. 


II. MASS-INDEPENDENT REACTION RATES 
AND DIFFUSION COEFFICIENTS 

Here we consider the model with mass-independent re¬ 
action rates and diffusion coefficients. Starting with ED], 
all subsequent work on aggregation processes with a lo¬ 
calized source was done in the framework of this simplest 
model. (Studies of driven systems with three-particle ag¬ 
gregation [21] and of aggregation-annihilation processes 
[22] with a localized source also used this framework.) 
Furthermore, to avoid the challenge of treating infinitely 
many coupled non-linear PDFs, the process was often in¬ 
vestigated only on the level of coalescence [SH] [23] when 
the distinction between individual clusters is ignored, so 
there is just one reaction rate and just one diffusion co¬ 
efficient. 

Even for aggregation processes with mass-independent 
reaction rates and diffusion coefficients some basic quan¬ 
tities haven’t been analyzed. For instance, the behav¬ 
ior of the monomer density was apparently unknown; we 
will show that it decays algebraically, ci ~ ^ 

far away from the source (in three dimensions). To de¬ 
rive this and other results we start with the full time- 
dependent equations 

— = V^Cfc + ^ acj - 2cfeC -f J6k,iS{r) (4) 

i+j=k 

Here c(r,f) = Cfc(r, 0 is the total cluster density. 

Equations Q follow from (|^ when Kig = 2 and Dk = 1. 
This choice makes formulas less cluttered, and it does not 
limit generality since it essentially amounts to choosing 
the units of length and time. Indeed, consider the general 
case of arbitrary mass-independent reaction rates and dif¬ 
fusion coefficients: Kig = 2K and Dk = D. If L denotes 
the unit of length and T denotes the unit of time, the 
dimensions of the parameters D and K are [D] = L'^jT 
and [K] = L'^/T. Thus we can use [K/ as a 
unit of length and j as a unit of time [24] . 
Measuring in these units is equivalent to setting K = 1 
and D = 1. 


To appreciate the relative simplicity of Eqs. Q we em¬ 
phasize two its consequences. First, we notice that the 
mass density M{r,t) — ^^^(r, t) satisfies a closed 

equation, the diffusion equation with a localized source: 

^=^^M + JSir) (5) 

Equation (§ is already valid when the diffusion coeffi¬ 
cients are mass-independent (the reaction terms do not 
affect the evolution as the aggregation process is mani¬ 
festly mass-conserving). Second, the total cluster density 
also obeys a closed PDE 

|^ = V2c-c2 + J5(r) (6) 

The validity of Eq. ^ relies on the mass-independence 
of both the reaction rates and diffusion coefficients. 

Equation is a non-linear PDE which has not been 
solved, but the asymptotic behaviors which are valid far 
away from the source are known |20j . We now outline 
these asymptotic behaviors and briefly discuss the match¬ 
ing with near-source behaviors. These results will be use¬ 
ful once we turn to the analysis of the full system Q . 

Let us start with the most relevant three-dimensional 
case. In the long time limit the total density should be¬ 
come time-independent. In the stationary regime in three 
dimensions we need to solve an ODE 


1 d / 2 dc 
dr V dr 


+ J5{v) = 0 


(7) 


The source term vanishes for all r > 0. Away from the 
origin, Eq. Q reads c" -I- ^c' = and admits a simple 
exact solution 


c = 


2 


( 8 ) 


This is just a special solution, not the general solution 
of the second order differential equation Q. The source 
term actually prescribes the boundary condition, 


lim 47rr^ ^ — J, 

r-fO dr 


(9) 


and (8) disagrees with this boundary condition. Equa¬ 
tion (8) yields the leading asymptotic far away from the 
source. Near the source, one anticipates that the reaction 
term is negligible and hence c = J/(47rr). One straight¬ 
forwardly finds further corrections: 


J 

47rr 



2 

ln( Jr) -I- . .. 


( 10 ) 


Equating the two asymptotic behaviors, r~^ ~ Jr~^, 
we estimate the location of the crossover region: r ^ J~^. 
The asymptotic ^ [resp. ( [l0| ] is therefore applicable 
when r J~^ [resp. r «C 
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The location of the crossover can be alternatively es¬ 
tablished without knowing the asymptotic behaviors. In¬ 
deed, the governing equation 0 and the boundary con¬ 
dition 0 are invariant under the one-parameter trans¬ 
formation group 

r —>■ Ar, c —A“^c, J —>■ J (11) 

Therefore the solution to Eq. Q must have the form 

c(r) = J2e(i?), R=Jr (12) 

The scaled density C(i?) is universal (that is, independent 
on J) and it satisfies 


1 d 

W Jr 


R^ 


de 

dR 


-e^ + (5(R) = 0 


(13) 


The crossover between near and far-from-source region 
must occur when i? ^ 1, i.e., when r ^ J~^. 

In terms of the scaled distance, the asymptotic behav¬ 
iors of the scaled density are 




l/(47ri?) 

2/i?2 


1 

i?> 1 


(14) 


The universality of the asymptotic behavior Q is 
remarkable—it does not depend on the strength of the 
source. Thus the details of input play no role at all 
as far as the leading asymptotic behavior of the total 
density is concerned; near the source the details become 
relevant. In the most general case when various clus¬ 
ter species are injected, so the source is characterized by 
{Jfc| k = 1, 2,...}, the far-away asymptotic is still given 
by the same formula Q, while near the source we can 
use ([^ with J = J2k>i Jk- 

Consider now our driven aggregation process Q in ar¬ 
bitrary spatial dimension d. The analog of ® reads 


c = 


2(4 - d) 


(15) 


The analytical solution is again unavailable. The invari¬ 
ance of the governing equation allows us to write 


c(r) = j2/(4-<i)e(i?), R = ji/(4-d)^ 

The asymptotic behaviors of the scaled density are 


(16) 


(l/[{d-2)n,R 

[2(4-d)/R2 


d-2l 


i?< 1 
i?> 1 


(17) 


where = 27r‘^/^/r(d/2) is the surface “area” of the 
unit sphere in d dimensions. 

These results are valid only when d < 4. The vanishing 
of the amplitude in (15) at d = 4 suggests a possible 


logarithmic correction, and it is indeed easy to establish 
m the asymptotic behavior far away from the source 


c ~ 


,2 


(18) 


Note that the cluster density is still independent on the 
strength J of the source. The crossover from the far away 
asymptotic ( [T^ to the near-source behavior c ~ 
occurs when r oc exp(87r^/J). 

For d > 4, the reaction term is negligible far away from 
the source ^D]. Close to the source, the density is found 
by solving V^c -I- J5{v) = 0 leading to 


c ~ 


{d-2)ndrd--^ 


when d > 4 


Far away from the source, the governing equation is 
asymptotically the same, but the source should be renor¬ 
malized. The solution thus has the same form (19), with 


a renormalized source strength Jeff- The reaction terms 
are important when r ~ 

An extra care is required in dimensions d < 2. The rate 
equations (|^ are valid above the critical dimension. For 
the constant-kernel aggregation, the critical dimension is 
dc = 2 (see e.g. [H]). In two dimensions, the correct 
behavior is obtained after a mild modification of the rate 
equations, e.g., ([^ should be replaced by 


F)c 

— = V^c — (?I ln(l/c) -I- J8(v) 
at 


from which 


81 nr 


c ~ 


( 20 ) 


( 21 ) 


V-i 

d= 1 

Inr 

d= 2 

^-2 

CO 

II 

r“^(lnr)“^ 

II 

jr-id-2) 

d > 4 


far away from the source. 

Collecting the asymptotic behaviors far away from the 
source and adding the one-dimensional asymptotic (see 
[20]), we get 


( 22 ) 


The total number of clusters IN'(t) = J drc(r,t) also 
exhibits amusing behaviors. The density quickly vanishes 
beyond the diffusive length r ~ Jt. Thus we use ( |2^ and 
perform the integration up to r ~ Ji to hnd |25j 

dn t d = 1 
(Int)^ d = 2 

N - <1 d = 3 (23) 

t/lnt d = 4 
Jt d > 4 


[22]) for the total cluster density were 
In Eqs. (22)-(23) we additionally 


Spatial behaviors 
established in Ref. 
emphasized a rather paradoxical lack of dependence on J 
in d < 4 dimensions. Thus the entire process is driven by 
the source, yet the leading behavior of the total number 
of clusters in the system is insensitive to the strength J 
as long as the source is non-vanishing. 

Little is known on individual cluster densities. We now 
establish some of their properties. 
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A. Monomer Density 


The asymptotic behaviors of the monomer density in¬ 
terestingly vary with spatial dimension: 


Cl 


V-5 d=l 

< j-(\/T7-3)/2^-(^-H1)/2 ^^3 

J-ir-2(lnr)-2 d = 4 

d>4 


(24) 


leading to the announced asymptotic Ci ~ r The 
dependence on J appearing in ci ^ j-(V 2 -i} jg 

established using the invariant form of the solution: 

c{r) = je{R), ci{r) = JCiiR), R = ry/j (31) 

Treating the dimensionality d as a continuous parame¬ 
ter we find that in the general case 2 < d < 4 the density 
of monomers far away from the source obeys 

Cl ~ (32) 


In one dimension, the analysis H OI HD requires separate 
techniques (outlined below in Sec IID) as the rate equa¬ 


tion framework is inapplicable. The higher-dimensional 
behavior is simple to appreciate. Indeed, above four di¬ 
mensions clusters do not “see” each other far away from 
the source and hence Ck ~ is valid for all k. The 

four-dimensional case is most tractable and we shall de¬ 


termine the entire cluster mass distribution (Sec. IIB|. 


In two and three dimensions, the exponent describing the 
large distance behavior is irrational. 

We now deduce these unusual behaviors from the gov¬ 
erning equations. When r \/i, the monomer density 
is stationary. Thus in three dimensions we must solve 


V^ci - 2cic -b J6{r) = 0 


(25) 


Recalling that c = 2/r^ in three dimensions far away 
from the source, we recast (251 into 


1 d 

r2 


dci 

dr 


4ci 

r .2 


(26) 


The general solution of this linear ODE is a combination 
of pure algebraic functions, with /3± = — 4 ± ^a/TT. 
To ensure that the density does not diverge as r —>■ oo, 
we pick Cl = Thus we arrive at the an¬ 

nounced irrational exponent. It is impossible to deduce 
the dependence of the amplitude A on the strength J of 
the source from a linear Eq. (26). However, we can de¬ 


termine the amplitude on the general grounds, viz. using 
the invariance of the governing equation (25). Equations 
0 and ( |25[ ) are invariant under the transformation group 

r —^ Ar, c —^ X ^c, Ci —^ A ^Ci, J —y A (^^) 

and therefore 

cir) = j'^e{R), ci(r) = j2ei(R), R = Jr (28) 

Erom Cl oc 7•-(^W7+l)/2 ^ j^-{VT7+i)/2 leading 

to Cl oc J-(^-3)/2, 

In two dimensions, we must solve 


V^ci = 


2cic 


ln(l/c) 

Using (21) and keeping only the leading term we get 
1 d 


(29) 


r dr 


dci 

dr 


8ci 


(30) 


The exponent /3(d) is found by solving V^ci = 2cic, or 
1 d / ^_i dci A 4(4 — d)ci 

j.d-1 dj- y dr I 


to yield 


/3(d) 


d - 2-b y(d - 2)2-b 16(4 - d) 
” - — • 


Another exponent 


a(d) = 


m - 2 

4-d 


(33) 


(34) 


is determined using the invariant form of the solution, 
ci(r) = j2/(4-^)ei(i?), R = 
together with Ci ~ 


B. Mass Distribution in d = 4 


The marginal case of d = 4 is the simplest and all den¬ 
sities c/c(r) can be determined analytically. Let us start 
with monomers. In the stationary regime, they satisfy 
V^Ci = 2cic, so in four dimensions 


1 d / 3 dci A 4ci 
r3 dr V dr / r'^t 


(35) 


where we used (21) and the shorthand notation £ = Inr. 
Solving (35) gives the leading behavior 


Cl = 


A 

r2 p. 


(36) 


with yet undetermined amplitude A. Similarly for dimers 
the governing equation reads 


1 d / 3 dC2 A _ 4 c2 

r3 dr V dr ) r'^i 


(37) 


Plugging (36) into (37) and solving for the density of 
dimers we find 


A 

^2 - 


(38) 


























5 


with the same amplitude as in (36). Solving for a few 


more densities we guess the structure of the solution: 


Cfe = ^ Ck{(.) 


Plugging this ansatz into 


J.3 


dck 

dr 


CiCj = 2cfeC (40) 


we arrive at 


^ = E C.C,-2CkC 


(41) 


i+j=k 


where C = J2k>i^k is asymptotically equal to £ ^. 


Equation (41) is identical to the solvable Smoluchowski 


for the constant-kernel aggregation—the variable £ = \nr 
plays the role of time. The solution reads 


„ A f Ak 

ft=pexp -- 


(42) 


-1 


and it manifestly satisfies the sum rule J2k>i = £ 

To fix A we recall that the mass density M = J2k>i ^^k 
satisfies (§ which admits an asymptotically stationary 
solution when d > 2. For d = 4, we have 


M = 


J 




(43) 


which in conjunction with (39) yields the sum rule 




k>l 




fixing the amplitude A = Stt^/J. 

Combining previous findings we write the mass distri¬ 
bution in four dimensions 


167r^ 1 


Cfc — 


J r2£2 


exp 


k 

IF ~£ 


(44) 


C. Mass Distribution in d = 3 

Let us immediately seek the mass distribution in the 
scaling form 


Cfc(r) = r ^ F(x), x = - 

r 


(45) 


In the scaling regime, clusters are typically large, so we 
must require r ^ J~^ to assure that the reaction term 
dominates over the source term. 


The choice of the scaling form (45) is suggested by the 
sum rules 


E ’ E (46) 

k>l k>l 


The second sum rules merely gives the mass density 
which is the solution of Eq. ([^ in d = 3. The sum rules 
(46) are manifestly obeyed and they imply two integral 


( 39 ) relations on the scaling function F{x)-. 


poo poo 

/ dxF{x) = 2, / dxxF{x) = 

Jo Jo 


J 

47r 


(47) 


Plugging (45) into the governing equations yields 

ad^F ^ dF 

-I- 6x- -1- 2F 

ax^ ax 


[ dyF{y)F{x-y) = 0 (48) 
do 


The integral is the convolution. Therefore we make the 
Laplace transform, $(s) = dx e~^^F{x), and recast 
(48) into an ordinary differential equation 


s2(j)" _ 25$' - 2$ -I- 4)^ = 0 


(49) 


Here the prime denotes the differentiation with respect 


to s. The sum rules (47) become 


$(0) = 2, $'(0) = -^ 


(50) 


It appears impossible to find an analytic solution of 


(49) subject to (50). An asymptotic analysis gives 

P ^ .r(vT7-5)/2 


0 






„-3 


167r2 


(51a) 

(51b) 


The cc —>■ 0 behavior of F{x) can be extracted from the 
s —>■ oo behavior of 4>(s). Since $(oo) = 0, we can drop 
the non-linear term in (49) when s ^ oo. Solving the 
linear equation we get 4) ^ s-i'/Tf-3)/2^ from which we 
deduce (51a). The cc — )■ oo behavior of F{x) is incapsu- 


lated in the s —)■ 0 behavior of 4>(s). The latter is deduced 


from (49) and (50) by straightforward manipulations: 


4? = 2 — -— s — 

47r 327r2 


In s - 


This small s expansion yields (51b). 


Combining the asymptotic behaviors (51a)-(51b) with 


( |45[ ) we arrive at the asymptotic behaviors of the mass 
distribution in three dimensions: 


Ck{r) ~ A:(^-5)/2 ^-(vT7+i)/ 2 ^ ^ (52a) 

k^ r (52b) 


Ckir) ~ k ^ 


167r2 

We emphasize that these results are valid only far away 
from the source: r J~^. 


D. Mass Distribution in d = 1 

A lot is known in one dimension [iniiiiiiiniiii]- Here 
we recall a few basic results and provide more details for 
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the behavior at the origin which is rather peculiar: The 
densities turn out to be are finite (in higher dimensions, 
they are infinite). 

In the stationary regime, the governing equation for 
the total cluster density 

+ J5{x) = 0 (53) 

is solvable [20] 


c(x) 


6 



(54) 


The structure of (54) agrees with the general form (^16 


and it gives the explicit form, C(X) = 6/[|X| + (24)^/^] 
with X = of the scaled density in one dimension. 

The stationary cluster densities are encapsulated in the 
generating function [iniiiniiii] 


- "'■)«<") = (M + f(i-.)--/3p <“> 


The densities remain finite at the origin. Specializing 
(55) to X = 0 and expanding in powers of z one finds 


' Uj r(i)r(fc + i) 


(56) 


The above results, e.g. the exact total density (54), or 
following from (55) exact cluster densities 


ci(x) = 


4i 

(N+^)^ 


C2{X) = - 




3(|x|+^)3 (|x|+£)4 


56 

Csixj = — 


8 £3 


27(|x|+€)3 3 (|x|+^) 4 9(|a;|+^)5 


etc. disagree with asymptotic predictions, e.g., with far 
from the source asymptotic behaviors c \x 


^ and 


Cl 3^ |x| 3 appearing in (22) and (24). This is because 


the above exact results have been established in the rate 


equation framework, while in Eqs. (22) and (24) we cited 


predictions for the truly one-dimensional process. For our 
aggregation process, the rate equations are valid only for 
d > 2. This does not imply, however, that the above 
exact results are useless. Indeed, if the aggregation pro¬ 
cess occurs in three dimensions and new particles are 
injected on the two-dimensional interface, the rate equa¬ 
tion framework is valid. 

We now illustrate the techniques which are used for the 
treatment of the truly one-dimensional process. These 
techniques are rather special, see mi Hanoi HD, and they 
are ultimately based on the assumption that the diffusion 
coefficients are mass-independent. Another crucial pos¬ 
tulate is the assumption that clusters are size-less (point 
particles if the process occurs on the line, or each cluster 


occupies a single lattice point). We also note that in one 
dimension, it is rather natural to examine particles un¬ 
dergoing biased motion, and in this situation one often 
needs even more complicated methods, see [211120]. We 
always consider symmetric motion. 

A powerful treatment of the one-dimensional aggrega¬ 
tion process with mass-independent diffusion coefficients 
is based m on functions Pk{x,y;t), where x < y and 
fc = 0,1,2,_ By definition, Pk{x,y]t) is the probabil¬ 

ity that the total mass contained in the interval [x^y] at 
time t is equal to k. Only the total mass is taken into 
account: For fc = 3, for instance, there could be three 
monomers, or one monomer and one dimer, but these 
details (and the precise locations of the clusters in the 
interval [x,?/]) do not enter into the description. The 
probabilities Pk{x,y,t) evolve according to 

pJ[Pk-i{x,y,t) - Pk{x,y]t)]e{y)6{-x) (57) 


with Heaviside step functions assuring that the source 
contributes only when x < 0 < y. Equations (57) are 


linear and solvable, although the solution is rather cum¬ 
bersome. In contrast to higher dimensions where the 
governing PDFs are non-linear and non-solvable, in one 
dimension one can obtain exact results for the full time- 
dependent behavior. 

Let us look at the stationary regime and additionally 
focus on a symmetric (in space) subset of probabilities: 
Qk{x) = Pk{—x,x]t = 00 ). Equations (57) give 


0 — 2 ——I- J[Qk-i — Qk] 


(58) 


The generating function Q{x,z) = J2k>o ^’^Qkix) satis¬ 
fies 

2^ = J{l-z)Q (59) 

Solving this equation subject to the boundary conditions 


Q{x = 00 , z) = 0, Q{x = 0, z) = I 


we obtain 


Q{x,z) = exp 



(60) 


(61) 


The densities at the origin are Cfe(O) = Hence 


^z^Cfc(O) 

k>0 


1 dQ 


2 dx 


x—0 



from which 


Cfe(O) 


r(fc-i) 

327 r r(fc-kl) 


(62) 


(63) 


Pushing these calculations one can determine the spatial 
distributions Ck{x). 
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III. BROWNIAN COAGULATION WITH 
SPATIALLY LOCALIZED SOURCE 

In this section we study Brownian coagulation driven 
by a localized source of monomers. We begin with the 
classical Brownian coagulation when clusters are spher¬ 
ical [Ml HU EH 12^ . In three dimensions the reaction 
rates are given by p|) and the diffusion coefficients are 


Dk = 


d-2 


, Dk ^ 


(64) 

-(d-2) 


Generally Kij ~ {Di -|- Dj){Ri + Rj) 
in d > 2 dimensions, so (ignoring numerical factors) the 
reaction rates are 

A,,, = (65) 

and the diffusion coefficients are 


Dk = 


( 66 ) 


for spherical clusters undergoing Brownian coagulation 
in d > 2 dimensions. 

We then briefly explore what can happen when clus¬ 
ters are fractal. In three dimensions we still can use the 
general prediction Kij ~ {Di + Dj){Ri + Rj) of the reac¬ 
tion rate theory OlllII] and the Stokes-Einstein relation 
Dk ~ 1/Rk- The characteristic size of a fractal cluster 
grows with its mass as Rk ^ k°‘, where a = 1/Df is 
the inverse fractal dimension of clusters. Thus in three 
dimensions, the generalized Brownian kernel is 

= (*“+J“)(^■“+r“) (67) 

and the diffusion coefficients are 

Dk = ( 68 ) 

A. Classical Brownian Coagulation 

Consider the classical Brownian coagulation and as¬ 
sume that the densities becomes stationary in the long 
time limit. The governing equations read 


Dk^ C/j; -f- ^ ^ ( KijCiCj Ck ^ — 6 (69) 

i+j=k j>l 

for k > 2. The density of monomers satisfies 


UiV^ci - Cl ^ KijCj + Jd{r) = 0 (70) 

i>i 

This infinite system of coupled non-linear ODEs with 


reaction kernel (65) and diffusion coefficients (66) is an¬ 


alytically intractable. Let us try to establish major fea¬ 
tures using heuristic arguments. 


Summing all Eqs. J69| and (70) we obtain 


^ DkCk -IY. + Jd(r) = 0 (71) 


Let s = s(r) be a typical mass. The typical diffusion 
coefficient is Dg = Invoking scaling we estimate 

J2k>i^kCk ~ The kernel (65) has homogene¬ 
ity index zero since it satisfies = Kij. There¬ 

fore we replace the Brownian ( |65[ ) by constant kernel in 
scaling estimates. This gives j>i ~ so 

Eq. ([tT]) leads to relation 


r-2^-l+2ldc^^2 


(72) 


To establish another relation we multiply (69) by k, 
sum over all fc > 2, and also add (70) to give 


Y kDkCk + J5 {y) = 0 


k>l 


from which 


Yl ^^kCk — (j _ 


j 


1 


fc>i 


(d- 2)D.a r 


d-2 


(73) 


(74) 


where = 27r‘^/^/r(d/2). On the level of scaling we es¬ 
timate kDkCk Therefore Eq. (74) implies 




J 

„d—2 


From relations (72) and (75) we deduce 

s - r‘^jd/{4-d} 




(75) 


(76) 


Varying d and looking at the dependence on J we realize 
that the above results are valid when 2 < d < 4. The 


lower bound is already evident from Eq. (74), the sta¬ 


tionary solution (74) exists only when d > 2. The upper 
bound emerges from (76). For d > 4 there is no scal¬ 
ing, e.g. a finite fraction of monomers escape to infinity 
without participating in aggregation events. 

The total number of clusters 3Nl(t) = J drc(r,t) is esti¬ 
mated by using (76) to give 


N(t)~ /’^drr‘^-ic(r) ~ ^ 

Jo Jo ^ 


The integral varies logarithmically with time. One can 
make a more precise estimate of the cutoff length R. 
We write R ~ '/Dt with D ^ g-i+ 2 /d (66)] 

leading to R^ ~ Further, Eq. (76) yields 

s ^ J‘^Z(4-d)_ Fi-om these two relations we obtain 


i?- 


j-{d-2)/{4-d} j 


1 1/d 


(77) 


Thus 


k>l 


^i>i 


N(t) ~ J-(d-2)/(4-d) ^^^j-(d-2)/i4-d) 

In three dimensions our main scaling results read 
c~J-V-3, s-j3r3, N-J-Mn(t/J) (79) 






















These scaling laws imply that the mass distribution in 
three dimensions has the scaling form 

Cfc(r) = x = (80) 

It appears impossible to extract quantitative informa¬ 
tion even about limiting behaviors of the scaled mass 
distribution F{x). For instance, one would like to estab¬ 
lish how the monomer density decays with distance. We 
estimate 


j>i f>i 


and hence Eq. (701 gives V^ci ^ r~^ci. This equation 
is invariant under the scale transformation, r —>■ Ar, and 
hence the solution is algebraically decaying: ci ~ r~^. 
To determine the decay exponent, however, we need to 
know the precise numerical factor which cannot be ob¬ 
tained in the framework of the heuristic approach used 
in this section. 


B. Generalized Brownian Coagulation 


If clusters are fractal, the proper reaction rates and dif¬ 


fusion coefficients are given by (67)-(68|. Various frac¬ 


tal dimensions, and respective values of the parameter 
a = l/Df in Eqs . ([^ "(68), appear in applications. 
In polymerization |31) if emerging polymers are linear 
chains, the value a = I describes the situation when the 
polymers are stiff (so their fractal dimension is Df = 1). 
Ideal polymer chains (also known as polymers in the 6 
solvent) are essentially random walks, so Df = 2 and 
a = |. Polymers in a good solvent are essentially self¬ 
avoiding walks, so a « |. 

Assuming that the stationary regime is reached, we 
perform the same analysis as in the previous subsection 
and obtain Eq. ( |7T| ) leading to ~ c^. Equations 

(73)-(74) still hold and from the latter we deduce another 
relation: ~ J/r in three dimensions. Thus the 

cluster density and the typical mass are 

j-a/(l-2a)^-(2-3a)/(l-2a)^ ^ ^ (j^)l/(l-2a) (gj) 


Using the former result we find that the total number of 
clusters behaves according to 




J-^ X 


{ (j2i)(l-3a)/(2-3a) 

ln(t/J) 

const 


a< 3 


3 < a 


(82) 


In three dimensions Df < 3 and hence a > so in 
the physically relevant cases we anticipate that the total 
number of clusters either grows logarithmically with time 
(a = |), or saturates (| < a). 

Equation (811 shows that the exponents diverge when 


from below. Eurther, the predictions Eq. (81) for 


a > ^ are dubious. These are the indications that the 
stationary regime is never reached when a > |. A sim¬ 
ilar phenomenon, namely that three-dimensional Brow¬ 
nian coagulation with reaction kernel (67) never reaches 
a steady state when a > |, has been established in the 
case when the source was homogeneous [TB]. Mathemat¬ 
ically, the evolving system with homogeneous source is 
much simpler as one deals with ODEs, but even in that 
case the numerical analysis is arduous [HI [19] as the non¬ 
stationary cluster mass distribution is very tricky—there 
is e.g. a boundary layer region describing the small mass 
tail of the cluster mass distribution. 

Thus we don’t understand the behavior of the gener¬ 
alized Brownian coagulation processes with a localized 
source in the region a > ^ which includes a = ^, « |, 1 
arising in the context of polymerization. Let us also think 
more carefully about the region | < o < | where clus¬ 
ters are fairly compact, 2 < Df < 3. In this region the 
total number of clusters remains finite. The rate equa¬ 
tion framework is deterministic, so it is applicable only 
to systems with infinitely many clusters. For finite sys¬ 
tems there are always fluctuations. If the total number 
of clusters diverges as t —>■ oo, the relative magnitude of 
fluctuations vanish (as I/-v/N in most cases). For the clas¬ 
sical Brownian coagulation fluctuations vanish, although 
extremely slowly, viz. as (lnt)“^/^; for the generalized 
Brownian coagulation processes with a localized source 
fluctuations never vanish in the region | < a < |. 


IV. SUMMARY 

Reaction-diffusion processes driven by localized input 
are often found in Nature and used in industry. Some of 
these processes which occur e.g. in electropolishing [55] . 
dissolution [51], corrosion [5S], and erosion [55] involve a 
few species of atoms. These processes are well understood 
as they are tractable mathematically [57114(1] . In other 
examples we have numerous interacting sub-species, e.g. 
clusters in aggregation [50] , aggregation-annihilation [55] , 
and mass exchange [IT] ; the analysis of these systems are 
much more challenging. 

Here we investigated aggregation processes driven by a 
localized source of monomers. When reaction rates and 
diffusion coefficients are mass-independent, a stationary 
regime is reached in the long time limit. In particular, 
we demonstrated that the densities of clusters of small 
mass exhibit an algebraic decay far away from 

the source. The decay exponent has irrational values in 
physically relevant dimensions: /3(3) = (vTf -I- I)/2 and 
/3(2) = v^. The understanding of this driven aggrega¬ 
tion process is quite complete in four dimensions where 
we presented the scaled mass distribution (again in the 
stationary regime). In three dimensions, the stationary 
scaled mass distribution satisfies a non-linear integro- 
differential equation. 

We also studied a Brownian coagulation process driven 
by a localized source of monomers. We argued heuristi- 












9 


cally that the stationary regime is reached, and we estab¬ 
lished chief scaling laws for a few major quantities, see 
Eq. (79), in the classical case when clusters are spherical. 
If clusters are polymers undergoing Brownian coagula¬ 
tion, one is naturally led to the reaction rates (67) and 
diffusion coefficients (68) with a > ^. In these situa¬ 
tions the stationary regime is apparently never reached. 
More precisely, we provided heuristic evidence that the 
assumption that a stationary regime is reached leads to 
unphysical behaviors. Our heuristic estimates relied on 
scaling assumptions, so the evidence is far from strong. 


In particular, the marginal case of a = ^ is known to 
be notoriously subtle for aggregating systems with a ho¬ 
mogeneous source [32], and this may be also true in our 
case. An obvious challenge is to explore the generalized 
Brownian coagulation processes in the region a > ^. We 
also showed that even when ^ < a < |, the deterministic 
rate equation framework provides limited insight for the 
generalized Brownian coagulation process driven by a lo¬ 
calized source. The reason is the finiteness of the total 
number of clusters implying the lack of self-averaging. 
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